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Explicit expressions for magnetoresistance R of planar and layered strongly in- 
homogeneous two-phase systems are obtained, using exact dual transformation, 
connecting effective conductivities of in-plane isotropic two-phase systems with 
and without magnetic field. These expressions allow to describe the magne- 
toresistance of various inhomogeneous media at arbitrary concentrations x and 
magnetic fields H. All expressions show large linear magnetoresistance effect 
with different dependencies on the phase concentrations. The corresponding 
plots of the x- and .ff-dependencies of R(x, H) are represented for various val- 
ues, respectively, of magnetic field and concentrations at some values of inhomo- 
geneity parameter. The obtained results show a remarkable similarity with the 
existing experimental data on linear magnetoresistance in silver chalcogenides 
Ag2+gSe. A possible physical explanation of this similarity is proposed. It is 
shown that the random, stripe type, structures of inhomogeneities are the most 
suitable for a fabrication of magnetic sensors and a storage of information at 
room temperatures. 



PACS: 75.70.Ak, 72.80.Ng, 72.80.Tm, 73.61.-r 
1. Introduction 

It was established recently that new materials, such as oxides and chalco- 
genides, have often unusual magneto-transport properties. For example, the 
magnetoresistance becomes very large (the so called colossal magnetoresistance) 
in manganites or grows approximately linearly with magnetic held up to 
very high fields in silver chalcogenides :2.. The large linear magnetoresistance 
(LLMR) takes place in thin films of Ag2+$Se,Te in a wide region of temper- 
atures, from low (~ IK) till room temperatures (~ 300K). At the moment 
there exist two approaches in a theoretical explanation of a linear behaviour 
of the magnetoresistance. The first, a quantum one, is proposed by Abrikosov 
and is based on the quantum theory of possible changes of spectrum properties 
of semimetals or narrow gap semiconductors 8]- It can be applied, for exam- 
ple, for an explanation of low temperature properties of Ag2+&Te ( 4 ). The 
second approach is pure classical and is based on the importance of the phase 
inhomogeneities, which take place in these materials on small (till nanometer) 
scales, for an existence of the LLMR effect at moderate temperatures .2 . This 
approach is applicable at moderate temperatures and for systems where a mean 
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size or some characteristic size of inhomogeneities l c ^> lo (here lo is a free path 
length of the charge carriers). In the framework of the second approach Parish 
and Littlewood have proposed the network model constructed from conducting 
discs and have shown by its numerical simulation that the LLMR appears in 
this model, when the parameters of a disc's impedance (in particular, a mobility 
fi) are random and have a continuous wide distribution with (/i) =0 5 , 6] . This 
model is similar to the usual wire network model [7|, but, in order to describe 
a dependence of the effective resistivity on magnetic field, it uses as building 
blocks the four terminal discs instead of wires. 

In the framework of the classical approach there is another possibility to 
describe magneto-transport properties of inhomogeneous planar (or layered, in- 
homogeneous in planes, but constant in the direction orthogonal to planes) sys- 
tems in perpendicular magnetic field. It is connected with an existence of the 
exact dual transformation, relating the effective conductivity a e (and the effec- 
tive resistivity p e = a^ 1 ) of planar inhomogeneous self-dual two-phase systems 
without and with magnetic field |SJ |S] . The existence of this transformation is 
a direct consequence of the exact Keller - Dykhne duality of two-dimensional 
systems [TtjUTT] . 

In this letter, using this exact transformation and the known expressions for 
<r e of three inhomogeneous models with different random structures from |12| . 
we will give the explicit approximate expressions for the effective resistivity p e 
of self-dual two-phase systems applicable at arbitrary values of phase concentra- 
tions and magnetic fields and in a wide region of partial conductivities. We will 
present also the x- and 7J-dependencies plots of the magnetoresistance R(x, H) 
at some characteristic values, respectively, of magnetic field H or phase con- 
centrations. These plots unambiguously show the existence of the large linear 
magnetoresistance effect in these classical systems. A comparison of results, ob- 
tained here analytically, with the known experimental data on the magnetoresis- 
tance behaviour in silver chalcogenides Ag2+gSe demonstrates their remarkable 
qualitative similarity. A physical explanation of this similarity is proposed. A 
possibility of an application of random, stripe type, inhomogeneities for a con- 
struction of magnetic sensors and magnetic read-write technologies is indicated. 

2. Effective resistivity in magnetic field 

The effective conductivity of two-phase isotropic systems in magnetic field 
has the following form 

& = °ik = VdStk + vttik, e>d(H) = <r d (-H), cr t (H) = -ot(-H), (1) 

here 5ik is the Kronccker symbol, is the unit antisymmetric tensor. The 
effective resistivity p e is defined by the inverse matrix 

p = Pik = P&kk + Pt<kk, Pd(H) = pd(-H), p 4 (H) = -p t (-H), (2) 

where 

Pd = <Jdl{p d + of), p t = -VtKo\ + Ot). (3) 
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The effective resistivity p e and, consequently, p ed , p e t (we assume here that 
Pid > 0) for self-dual systems must be a symmetric function of pairs of partial 
arguments (f>i,Xi) and a homogeneous (a degree 1) function of pdi.u- For this 
reason it is invariant under permutation of pairs of partial parameters 

Pe(pi,X 1 \p 2 ,X 2 ) = Pe(p2,X 2 \pi,Xi). (4) 

The effective resistivity of inhomogeneous systems must also reduce to some 
partial pi, when Xi = 1 (i = 1, 2). 

In our previous paper we have obtained explicit expressions for effective 
conductivities of planar inhomogeneous self-dual systems in magnetic field and 
have shown that they have properties qualitatively compatible with the large 
linear magnetoresistance effect. Since the people usually mesure in experiments 
the effective resistivities and such value as the magnetoresistance, we consider in 
this letter the properties of the p e d{x, H) and the magnetoresistance R(x, H) of 
strongly inhomogeneous two-phase planar (and layered) systems. The investiga- 
tion of properties of the classical Hall effect in these systems will be done in the 
subsequent paper. The magnetoresistance of two-dimensional inhomogeneous 
media is determined as follows 

R(x, H) = (p ed (x, H) - p ed (x, 0))/p ed (x, 0). (5) 

One needs to note here that there is another definition of the magnetic resistance 
in the literature, which differs from (5) by the normalization 

Ri{x,H) = (p ed (x,H)~ p ed (x,0))/p ed {x,H). (5') 

We will use here a definition (5), since it allows to show in a more obvious 
way the LLMR effect. For isotropic systems without magnetic field p ed (x,Q) = 
Peo = (feo)" 1 , where a e o is the effective conductivity at H = 0, and, for self- 
dual systems, p e Q has the same functional form as a e Q. Then one can write R 
as 

R{x,H) = a e0 p ed (x,H) - 1 = g e0 2 e 2 - 1. (6) 

a ed + a et 

It follows from (5) that R(x, 0) = and can be nonzero only for H ^ 0. Thus, 
it describes only transport properties connected with magnetic field. Note 
that for large values R can be approximated by <r e0 p ed . In order to describe 
the dependence of the R(x, H) of inhomogeneous system on magnetic field 
one must know also the functional dependences on H of partial conductivities 
(Tid(H-), (Tjt(H), (i — 1,2). They usually can be approximated by the standart 
(metallic type) formulae [TT] 

<7 ld (U) = ^ CTft(H) = ^^2 ' 0i=Vi H > i = h 2 , (7) 

where cr^o are the partial conductivities of phases without magnetic field, pi = 
eiTi/rrii are the corresponding mobilities, rj,e,,mi are, respectively, carrier's 
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times of life, charges and masses. Thus R(x, H) depends also on 4 partial pa- 
rameters: crj, Since a e d,t are homogeneous functions of cr^o and the mobilities 
enter always in a combination with magnetic field, really R(x, H) depends (be- 
sides x and renormalized magnetic field H' = piH, below we will assume that 
fi% = 1), on two inhomogeneity parameters: z, connected with an inhomogeneity 
of conductivities, and rj, connected with an inhomogeneity of mobilities. They 
may be chosen as the corresponding ratios (we choose also that 020 /cio < 1) 

2 = 0-20/0-10 (0 < z < 1), J7 = /i 2 /Mi (-oo<7/<oo). (8) 

It is useful also to note that p d for homogeneous systems with conductivities 
from (7) has the form pid{H) = a^ 1 , i = 1,2, and, consequently, for these 
homogeneous systems 

R(x, H) = <r ed {x, 0) Ped (x, H) - 1 ► 0. (9) 

x— 0,1 

The equation (9) defines the boundary values of R(x, H) at x = 0, 1 for inho- 
mogeneous systems, satisfying the representation (7). 

3. Magnetoresistance of planar inhomogeneous systems. 

Recently, using exact "magnetic" duality transformation, connecting the ef- 
fective conductivities of inhomogeneous systems with and without magnetic 
field |3, we have obtained three explicit approximate expressions for the effec- 
tive conductivity a e of planar two-phase strongly inhomogeneous systems with 
different structures of inhomogeneities in magnetic field Two of them de- 
scribe systems with various real "bulk" inhomogeneities (i.e. having real 2D 
bulk inhomogeneities): the compact inclusions of regions of different sizes and 
forms ("random droplets") of one phase into another and the "random par- 
quet" structure constructed from square plaquettes with randomly distributed 
stripes of two phases ^31- The third one has the effective conductivity, ob- 
tained by a dual "magnetic" transformation from the known effective medium 
approximation (EMA) formula, based on the wire-network representation of 
inhomogeneous systems |Hj (we will name it the effective medium (EM) model) . 

The corresponding effective conductivities have the following forms 

a ed {{a}, {x}) = e f v ' a et ({a}, {x}) = c ) f ; 2 , (10) 

where a' d is the effective conductivity of the models without magnetic field, but 
it depends on transformed partial arguments a' id = (Jid/cai with a ai = <j it + a. 
The parameters of the "magnetic" transformation a, b' = b/c,c depend on the 
partial conductivities and have the following form 

\a 2 \ 2 -\a^±VB fa] 3 - \v a \*±VB 

a ± = 7T, v ) b ± = TT, \ ) c = -«> 

2{(T lt -(T2t) 2(au-cr2t) 

B = [(oit - <J2t) 2 + [aid - °2d) 2 }[(<Tit - <J2t) 2 + (aid + V2d) 2 }, (11) 
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where \ai\ 2 — a 2 d + af t , and, evidently, B > 0. The effective conductivities 
fed({c}) {%}) of the three inhomogeneous models without magnetic field have 
simple forms: 

(1) a "random droplets" model 

a ed ({a},{x}) = (a 1 yi(a 2 r*; (12) 

(2) a "random parquet" model 

^(M,W) = (^) 1/2 ( ^^ ) 1/2 ; d3) 

(3) an effective medium model 



<Ted(W}, {x}) = (X- i)(<7i - CT 2 ) + J (x - \) 2 {<Jl ~ ^2? + <Tl<T 2 . (14) 

All formulas (10,12-14) correctly reproduce boundary values of the effective 
conductivities as well as their exact values at equal phase concentrations x = 1/2 

P2J 



a ed {x = 1/2) = y/aid<T2dA, A 



1 



&lt — &2t 



a et {x = 1/2) 



&ld + &2d 
cr 2 t&ld + G\t<J 2 d 



1/2 



&ld + &2d 

This permits us to write the exact expression for R at x = 1/2 



R(l/2,H) = a ed (l/2,0) 



<Ted(l/2,H) 



a 2 ed (l/2,H)+a 2 et (l/2,H) 



(15) 



where 



a ed (l/2,H) 



V^IO^O 



A(H), 



A(H) 



l + H 2 



a et {l/2,H) = 



^(1 + H 2 )(1 + V 2 H 2 ) 

(l + ^H^ + zil + H 2 ) 
a 20 (l + ri)H 



(1 + rf-H 2 ) + z(l + H 2 



(16) 



One can obtain an asymptotic behaviour of R at high H, r/H S> 1. Since in this 
limit 



a ed (l/2,H) = 



\/&10& 2 Q 

H 



rj — z 



J] 2 + z 



, a et (l/2,H) = 



^20(1 +V) 

(rf + z)H' 
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then 



R(l/2,H) = cH + o(H), c= J -^-0-±±^ + o{H). (16') 

Thus, the exact value R(1/2,H) has a linear dependence on H at high H with 
the coefficient c, which becomes for homogeneous media as well as for inho- 
mogeneous systems with rj = z. It is interesting that c = 1 at r\ = — 1. Just in 
this case a e t(x = 1/2) = 0. At T] = 1 c reduces to a simple formula c = ■ 

Substituting (11) and (12-14) into (10) and then into (6), one can obtain 
the explicit expression for the magnetoresistance R(x, H) of the corresponding 
systems at arbitrary concentrations and magnetic fields. Since these formulae 
are rather complicate we will analyze their behaviour by constructing the plots 
of their dependencies on x and H. 

Using explicit formulae (6), (7), (10-14) we have constructed the plots of H- 
and cc-dependencies of R{x, H) for inhomogeneous systems, whose inhomogene- 
ity structures are similar to three considered models, for different values of 
the inhomogeneity parameter z and r\. The corresponding plots of H- and x- 
dependencies at z = 10~ 2 and rj = 1 are represented, respectively, in Fig.l and 
Fig.2 ((a)-(c)). 

One can see from fig.l and fig.2 that the behaviour of the magnetoresistance, 
though different for various expressions (i.e. various models) and saturated in 
some regions of x, has two common characteristic features: 

1) absolute values of R are very large, especially at peaks; 

2) all values in the regions with large values of R increase with a growth of 
H at relatively high H approximately linearly. 

Both these properties are the consequencies of the high H behaviour of the 
effective conductivities ^21 

fed ~ cr ei ~ l/H, (17) 

and are qualitatively consistent with the experimental results from |2l H*4l IT5] . 
which show a large magnetoresistivity and its approximately linear growth with 
an increase of H. Practically the same plots for R(x, H) are obtained for the 
case, when two partial phases have carriers with opposite charges, i.e. at r] = — 1. 
It follows also from (6), since a e t, sensitive to a sign of rj, enters in R only as 
squared. Some dependence on a sign of 7/ appears for the Hall resistivity (see 
next paper). 

More detailed analysis shows a number of interesting specific properties. Let 
us start with the "random droplets" model. In this case R has evident one peak 
asymmetrical structure positioned for the inhomogeneity parameter z = 10~ 2 in 
the region of concentration x ~ 0.37 with a maximum growing almost linearly 
and slightly shifting to the higher values of concentration and with a width 
narrowing with an increase of H. Note that the maximal value of R is much 
larger then the exact value (15). It appears that a decline of linear dependence 
also changes becoming maximal ~ 2.7 at the peak and H = 300. For higher 



G 




H 



(c) 

Fig.l. The plots of the //"-dependence of the magnetoresistance R(x,H) for three 
explicit expressions obtained above (respectively, (a),(b),(c)) at the inhomogeneity 
parameters z = 0.01, and r\ = 1 at some characteristic values of concentration x. 

values of inhomogeneity parameter z the peak slightly shifts to the lower values 
of x (for z = 10~ 3 it is situated at x w 0.3). 

Very remarkable pictures are obtained for the "random parquet" model. The 
magnetoresistance shows approximately constant large values, almost coinciding 
with the exact value at x = 1/2 in a wide region of concentrations and very sharp 
changes in small regions x, 1 — x <C 1 near the boundary concentrations x = 0, 1. 
In this model the exact value is a maximal one. The width of the plateau 
region, which is approximately symmetrical, slightly grows with a growth of H. 
For the concentrations in the plateau region R increase linearly with H with 
an approximately constant coefficient, which coincides with its exact value at 
x = 1/2 from (16). The width of the narrow regions near boundaries decreases 
with an increase of H and a decrease of z. 

For the effective medium model we have obtained very unusual pictures. 
R(x, H) shows strongly asymmetrical very narrow and high peak, situated just 
below x = 1/2, with a sharp decrease to negligably small values in the narrow 
region near x = 1/2 and long small tale in the region x < 1/2. The exact value 
of R is significantly smaller than the maximal one. 

In order to compare R(x, H) for different models in more detail we represent 
also the plots of their ^-dependencies at H = 50 in one Fig. 2(d). One can see 
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(c) (d) 

Fig. 2. The plots of the x dependence of the magnetoresistance R(x,H) for three 
explicit expressions obtained above (respectively, (a),(b),(c)) at the inhomogeneity 
parameter z — 0.01, and at the four different (dimensionless) values of magnetic field 
H : 1) 50, 2) 100, 3) 300 and 4) 900 (the corresponding plots go from the lower to the 
upper ones; on (c) H = 25, 50, 100, since in this case for higher H the peaks are very 
narrow and almost coincide), (d) shows a comparison of all R(x,H) at H = 100. 

from Fig. 2(d) the relative widths and maximal values of R for three models. The 
highest values R achieves for the "random droplet" and effective medium models, 
meanwhile the "random parquet" model has the widest region of relatively large 
values. At equal phase concentrations x = 1/2 all models give exact value 
(16,16'). 

4. Discussion and comparison with the experiments 

The obtained formulae and the constructed plots unambiguously show an 
existence of the large magnetoresistance effect with almost linear growth at 
relatively large magnetic fields in various planar inhomogeneous models. In 
different models it has different dependences on concentration. The analysis 
of these dependences and the structures of the considered models allows to 
conclude that the different inhomogeneous fluctuations have different influence 
on the magnetoresistance. The random compact (droplet like) fluctuations have 
the most strong influence in the range of concentrations 0.3 < x < 0.5. In this 
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region the randomly distributed droplets of higher conducting component serve 
as traps for charge carriers. They work the most effectively for concentrations 
in the region 0.3 < 0.4 where all droplets have finite sizes and rather rare 
intersections. For larger concentrations the large droplets appear. Though the 
carriers are confine inside them for a long time, they already can propagate on 
significant distancies due to their large size. The last effect induces a decrease 
of R. The almost linear growth of R with H means that the stronger magnetic 
fields help to confine carriers inside droplets. 

At the same time, the stripe-like inhomogeneities, as it follows from the re- 
sults for the "random parquet" model, begin to work as traps already at very 
small concentrations x, when they have a form of rare, randomly distributed, 
stripes of highly conducting component. If a length of these stripes is con- 
strained and they do not form any effectively large or dense cluster (what takes 
place in this model due to a construction), these stripes gives approximately the 
same contribution to R, very weakly depending on their concentration. Only 
in very narrow regions near x = and x = 1, where a number of highly or 
badly conducting stripes is very small, the rare randomly distributed stripes 
give sharp changes of R, which quickly saturate. We suppose that a similar 
behaviour will take place in the models with other randomly distributed stripe- 
like inhomogeneities. For example, we expect that such behaviour (maybe, a 
more smoother than in the original "random parquet" model) will take place in 
models constructed from plaquettes of different forms and sizes with isotropic 
discrete (or even continuous) distributions of stripe orientations. 

The most unusual behaviour shows the " magnetically" transformed effective 
medium model. It does not influence on R almost for all concentrations, except 
very narrow region near x = 1/2 and long tail with non large values in the region 
< x < 0.4. It can be explained by the tight connection of this model with the 
wire- type networks, which cannot be used for a description of magneto-transport 
effects, connected with bulk properties. Only in the narrow region near x = 1/2, 
where the large dense clusters with many "dead" ends appear, they give strong 
contribution to R. 

Now we can compare the considered models with the known results for silver 
chalcogenides (Ag 2 +sSe and Ag2+sTe), which show large linear magnetoresis- 
tance effect (see, for example, for different concentrations of Ag ^5]. As 

it was noted in silver chalcogenides have a tendency to a formation of differ- 
ent phases on small scales (including nanometre one). The additional Ag ions 
strive to situate on various defects and grain boundaries ^HIE]- The latter pre- 
sumably have a form of intersecting straight lines. When concentration of the 
excess Ag ions is small, they form randomly distributed stripes of different size 
along these boundaries |15|. This structure is very similar to the corresponding 
structure of the "random parquet" model. The dependences of p e Q and R(x, H) 
on x in the "random parquet" model at small x (the ^-dependences of p e g in all 
three considered above models can be found, for example, in |13| . where they 
are represented for corresponding cr e o , one only needs to change o\ , o-i into p2 
and pi respectively) are also very similar to the behaviour of p e o and R(x, H) 
observed in experiments 115| , where very small excess of Ag ions gives sharp 
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increase of R. Moreover, for the larger concentrations, in some finite region, 
the experiment shows an approximately constant behaviour of p e o and R, which 
changes again on, respectively, decreasing and increasing behaviour |15j . We can 
explain this behaviour as a result of a crossover from the presumably striped 
structure of excess Ag ions to the mixed structure, including also their compact 
islands ^H]. The latter forms a structure, which is similar to that of the "ran- 
dom droplets" model. The corresponding contribution into R is firstly small 
relative to "random parquet" model. But, as follows from our Fig. 2(d) (its en- 
velope in the region < x < 0.5 qualitatively coincides with Fig. 3 from ITS]), 
at approximately x c = 0.25 (the exact values of the crossover x c can depend 
on magnetic field and inhomogeneity parameters) a contribution of the "ran- 
dom droplets" model becomes larger than that of the " random parquet" model 
and again appears the increasing behaviour of R (or a decreasing behaviour of 
peo). It means that the structure of inhomogeneities with the excess Ag ions 
in Ag2+sSe(Te) is very similar to the structures of our models (at least, in the 
region < x < 0.5) and the latter can give the necessary formulae for a descrip- 
tion of the magneto-transport properties of silver chalcogenides in this region 
of concentrations. Moreover, our formulae give also numerical values, which are 
in a good agreement with the experimental ones (see our forthcoming papers). 

At the end we would like to note a possible importance of the obtained results 
for practical applications. As is known, a behaviour of R such as it takes place in 
the "random" parquet model and models related with it is the most suitable for 
magnetic sensors and information storage technologies (see, for example, |17|1. 
For this reason, the "random parquet" type models with random stripe-like 
structures of inhomogeneities show us what kind of inhomogeneous structure 
one needs to fabricate sensors and other hightechnology devices, having small 
sizes and working at room temperatures. 

5. Conclusion 

Using three explicit approximate expressions for the effective conductivity 
of 2D isotropic two-phase systems in a magnetic field, we have obtained explicit 
expressions for the magnetoresistance R(x, H) of planar isotropic two-phase 
systems. The plots of the x- and ii-dependencies of R(x, H) at the different 
values of the inhomogeneity parameters z and rj are constructed. They show 
the evident large, almost linear, magnetoresistance effect. The behaviour of the 
constructed plots is very similar and qualitatively compatible with the experi- 
mental data for silver chalcogenides Ag2+sSe from [21 E3- A possible physical 
explanation of such behaviour is proposed. It is noted that the random, stripe 
type structures of inhomogeneities, are the most suitable for a fabrication of 
magnetic sensors and a storage of information at room temperatures. We hope 
that the constructed models and their modifications as well as the obtained 
results can be applied also for a description of magneto-transport properties of 
various real heterophase systems (regular and nonregular as well as random), 
satisfying the symmetries and having the similar structures, in a wide range of 
inhomogeneity parameters and at arbitrary concentrations and magnetic fields. 
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